subroutine comput_flux2
use var
implicit none

    call comput_Qy_m
    call comput_Qy_p
    call L_F_local_y
    !integer m,j,i
    !real(8) f(0:gridsx,0:gridsy,4)
    !real(8) local_max
    !real(8) f_m,f_p
    !
    !call comput_Qy_m
    !call comput_Qy_p
    !call L_F_local_y(Qy_m,Qy_p,f)
    !do m=1,4
    !    do j=0,gridsy
    !        do i=1,gridsx
    !            flux2(i,j,m)=f(i,j,m)
    !        enddo                
    !    enddo  
    !enddo
    
end subroutine